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Abstract: 

If string theory controls physics at the string scale, the dynamics of the early universe 
before the GUT era will be governed by the low-energy string equations of motion. Study- 
ing these equations, we find that depending on the initial conditions when the stringy era 
starts, and on the time when this era ends, there are a variety of qualitatively distinct 
types of evolution. Among these is the possibility that the universe underwent a period of 
infiation. A by-product of this analysis is the observation that it is often possible to erase 
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any evidence of a dilaton at late times. 
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In this paper we discuss the solutions of the lowest order string beta-function equations 
that represent homogeneous isotropic four- dimensional spacetimes. The basic ingredient of 
this calculation is therefore the massless states of the string, and for simplicity, we consider 
only the closed bosonic string. It is quite straightforward to extend our technique to other 
types of string theory. This type of solution is of importance for a number of reasons. The 
first is that it illustrates some of the possible solutions to the string equations of motion 
as a problem in its own right. One might think that the only solution to the bosonic 
string equations that is consistent with having spatial sections of spacetime being surfaces 
of constant curvature would be the well known example of 26-dimensional Minkowski 
spacetime. However, it is possible to trade off some of these extra dimensions in return for 
spacetime curvature, and this is precisely what we do here. 

Another reason for studying the problem is physically motivated. If string theory 
really controls physics at the string scale (presumed to be roughly Planckian although it 
could be at a significally lower scale), then from this era down to whenever the stringy 
symmetries are broken to yield the physics of (presumably) the GUT era universe, it seems 
reasonable to suppose that the dynamics of the universe will be governed by the low-energy 
string equations of motion. These are the equations to be studied in the remainder of this 
paper. 

Similar, but much more restricted results have previously been obtained by Antoniadis, 
Bachas, Ellis and Nanopoulos [1], Love and Bailin [2], Campbell, Kaloper and Olive [3], 
Tseytlin and Vafa [4] , and Tseytlin [5] . A number of other related references can be found 
in [5] . Their solutions are all special cases of our results which in some sense are a complete 
set of solutions that are consistent with symmetries of the spacetime that we are imposing. 

In the closed bosonic string theory, the long range massless fields are the dilaton the 
axion field strength Habc = [abc] = 69[ai?{,c] ? which is derived from the two- form potential 
Bab, and the graviton, or equivalently, the spacetime metric tensor gab- Condensates of 
these fields can be treated in a way consistent with the symmetries of the string and then 
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obey the j3 -function equations [6] . If we work in the string frame, these equations can all 
be derived from the action 

1 = dW'^e^ [c-R- {V<pf + ^HabcH'^'^^ (1) 

where R is the Ricci scalar of the metric gab, c is a constant, related to the central charge 
of the string theory, and the integral is taken over all of spacetime M. The scale of string 
physics is determined by by a', the inverse string tension. Provided that the spacetime 
curvature is small on the string scale, then equation (1) is a complete description of the 
massless modes of the string. However, for curvatures large on the scale of string physics, 
this action needs to be modified by terms in a'; however in the era we are considering, 
they are negligible. 

Variation of this action with respect to the metric yields the string analog of Einstein's 
equation 

Rab = ^a^b(f>+\HacdH^'' . (2) 

Variation with respect to Bab gives the axionic analog of Maxwell's equation 

VaH''^'' + Va(t>H''^'' = . (3) 

Finally, variation of the dilaton gives rise to 

c=R-{W(l)f-2ncj>+^HabcH''^' . (4) 

The constant c is given by ^^J^ for the bosonic string, and ^^^3^7"*^^ for the heterotic or 
superstring, where D is the dimensionality of spacetime and a' is the inverse string tension. 
However, c can be changed from these values by coupling some conformal field theory to 
the string, and so c can have any value in practice. One often studies critical string theory 
where c = 0. The reason for such a choice is that it enables Minkowski spacetime to be a 
stable ground state for the string. However, we are interested in cosmological solutions in 
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string theory, and so this is no longer a meaningful restriction. We regard c as an arbitrary- 
constant in what follows. 

The /3- function equations (1-4) describe physics as seen from the viewpoint of the 
string. However, they are not convenient for understanding gravitational phenomena, 
because the coefficient of the Ricci scalar depends on the dilaton field. An easy way to 
look at spacetime physics is to perform a conformal transformation on the metric so as 
to eliminate the dilaton-dependent term. This conformal frame is usually referred to as 
the Einstein frame, whereas the original one is termed the string frame. Such a conformal 
transformation yields a new metric gab given by 

9ab = e~'^gab (5) 

so that the action becomes 

/ = J dW'^ (ce-^ -R- \{Vct^f + ^e'^HabcH^'^^ . (6) 

We see from the new action that in the Einstein frame, gravitation is described by the 
minimal gravitational action but matter fields couple to gravity via the dilaton with various 
conformal weights. This is reminiscent of the situation in Brans-Dicke type of theories. 
Variation of (6) leads to the field equations in the Einstein frame. The equation of 
motion is 

□ = ce-'^ - ^e^-^ffabeiy"'" . (7) 
The H equations of motion are 

V[ai?6cd]=0 , (8a) 

VaH""^^ - 2Va0ii'"'" = , {8b) 

and the Einstein equation is 

Rab - IrRgab = ^ (^a(l>^b4> " l:9abi'^4>y 



+ [HacdHf - IgabHcdeH'^"'^] ■ (9) 

+ ^cgabe~^ 



Thus the Einstein tensor couples to the minimal energy-momentum tensor of the 
dilaton, and the minimal energy-momentum tensor of the H-field weig hted by e^*^. The 
central charge term c appears in a way reminiscent of the cosmological constant, but 
again with a weight factor now of e"*^. One might worry about the consistency of this 
set of equations. On general grounds, one knows that (4) follows from (2) and (3) as a 
consequence of the Bianchi identities. In a similar way, one can derive (7) as a first integral 
from (8) and (9). 

As is apparent from (9), the right-hand side looks like a conventional gravitational 
theory, and therefore one should not be surprised to discover singularities in the solutions of 
(9) where the curvature blows up. Such singularities are however, unlike general relativity, 
not necessarily unphysical. As far as string physics is concerned, one only needs to ask if 
quantities in the string frame blow up, since only then will the string be badly behaved. 
In other words, the types of singularity predicted by the singularity theories in general 
relativity do not necessarily cause breakdown of physics in string theory. 

The Universe on very large scales looks like a four-dimensional Friedman-Robertson- 
Walker spacetime with a value of k which cannot be observationally determined. Therefore 
on large scales the spacetime metric is 

ds^ = -dt'^ + a'^{t)dal (10) 

where dcr^ is the line element on a unit three-sphere, Euclidean 3-space, or the unit hy- 
perboloid depending on whether k = 1, 0, or —1 respectively, and a{t) is the cosmological 
scale factor. Consistent with the assumption of homogeneity and isotropy implicit in this 
metric, we choose the dilaton field to be independent of the special coordinates, and the 
axion field Habc = fif)^abc where Cabc is the volume form on the surfaces of t = const. 
Then the requirement that Habc be closed immediately implies that f(t) is a constant, /. 
We can then reduce the beta-function equations down to a set of three coupled ordinary 
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differential equations, 

0=-3- = 6 - +^--^+ 12 

where ' = d/dt. We now aim to characterize aU the solutions of (11) and (12). 

Although (11) and (12) are the equations in the string frame, it is straightforward to 
translate the solutions into the Einstein frame explicitly. In the Einstein frame, the metric 
is 

ds^ = e-'^ds^ = -dr"" + 6^ {r)da^ (13) 

where 

T = J dte-^'l' (14a) 

and 

6(r) = e-^^(^)a(T) . (146) 

To simplify the discussion of these equations, firstly we consider the special case of 
the critical string (c = 0) and a flat universe {k = 0). We can now eliminate / and find 
that 



0= = -12 - - 15-(/)-3(/)^ . (15) 

a \a/ a 

Perhaps the easiest way to find all the solutions to these equations is to change variables 
into 

V = 4> 

(16) 

X = o/o 

yielding 

f) = —12x2 — 15px — 3p2 

(17) 

X = 3x + 5px + P 

There are several simple analytic solutions that can easily [1,5] be found from (17). 



They are 



V- ^ X- ^ (18) 
t — to t — to 
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with any of 

i)A = B = 

u)A = -2/3 B = 1/3 

(19) 

iii) yl=l-v^ 5 = 1/^3 

iv) A = l + V3 B = -lV3 

Case (i) is flat Minkowski spacetime (in either the Einstein frame or the string frame) was 
first discussed by Antoniadis, Bachas, EUis and Nanopoulos, [1]. Case (ii) is a{t) = aot^^^ 
and (j){t) = (/»o — 2/2,lnt (setting — 0). If we transform this into the Einstein frame 
(00 = 0), then 

r=^t^/3 and 6 = feor'/^ • (20) 

In other words, case (ii) is identical to a radiation-dominated universe. It is singular at 
t = 0, even in the string frame, and unphysical since one finds that < 0. Cases (iii) and 
(iv), were originally found by Tseytlin, [5], have = and similar power law expansions 
and contractions. Case (iii) is a(r) = aot^^^ while case (iv) is a(r) = aQt~^/^. 

The remaining solutions can be explored by examining the phase-plane portrait using 
the fact that Eq. (17) results in 

dp ^ 12x' + Ihpx + 3p^ 
dx 3x^ + 5px + 

It is shown in figure (1), the region of physical solutions is bounded by the exact solutions 
of case (iii) and (iv) which are given by 

P=(-3tV3)x . (22) 

The evolution is in the direction of the arrows indicated. With the exception of the exact 
solutions (iii) and (iv), the spacetimes are nonsingular. There are two types of evolution 
depending on which side of the phase plane the trajectories begin. In both cases tracing 
the solutions backwards in time reveals that at earlier times they were repelled from the 
contracting version of either the exact solution (iii) or (iv). For solutions that expand 



rapidly at late times, the exact solution (iv) is an attractor. These type of solutions 
occupy the left-hand half of figure (1). The solutions that occupy the right-hand half of 
figure (1) are attracted to solution (iii) and on it to the (0, 0) point, i.e the universe comes 
to a halt asymptotically. In Figure (2) we plot a{t) and x(^) for two typical cases. 

We turn now to a more complicated situation in which the spatial sectors of the 
universe are not flat. In order to get a phase- plane picture similar to the one we got in the 
flat case , we change variables from t to 6 such that 



d9 _ 1 
dt a 



(23) 



Using ^ as a dynamical variable we now deflne x P to be: 

X=- P = <l>' , (24) 
a 

where ' = d/d9. The analog of Equation (17) can then be rewritten as 

P' = -UxP - 12x^ - 12A; - 3P^ (25) 

To check consistency, the solutions of Eqs. (25) must satisfy Eq. (11), which in terms of 
the new variables takes the form 

6xP + 6x^ + P^ + Qk = 2f/a^ . (26) 

We examine now the phase- plane properties derived from Eq. (25). Firstly we construct 
dP/dx and find that 

dP „ XP 

= _3 _| ± (27) 

dx 5xP + 4x2 4/c -t- P2 • V J 

The analytic solutions which bound the regions of the physical solutions are given by the 
requirement that = 0. From Eq. (26) it follows that 

6xP+Qx^ + P^ + 6k = . (28) 
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Therefore, in the asymptotic regime where k, the solutions which bound the physical 
solutions are given as before by Eq. (22) for the new P and %. However, the exact solution 
to Eq. (28) is 

P=-3x±V3x'-6A; • (29) 

Those, there is radically different behavior in the case k = 1 compared to k = —1 when 
one gets close to the origin of the x-P plane. We denote by (ASiii) and (ASiv) those lines 
that are asymptotic to solution (iii) and (iv), and turn now to examine the properties of 
the trajectories in a closed universe. 

For k = 1, physically valid solutions penetrate into the region that was forbidden in 
the case that k = 0. As in the flat case, there are two types of trajectories depending on 
which part of the phase-plane shown in Fig. 3 they occupy. All solutions are asymptotic 
to the lines given by Eq. (28). However, all the trajectories end up on the expanding part 
of (ASiv). These on the upper part of the (x, P) plane are repelled from the upper part 
of (ASiv), and are attracted to the collapsing part of (ASiii). The collapse slows down 
as the trajectories approach solution (ASiii). The trajectories route depend on whether 
they start below or above solution (iv). Trajectories that start from below go through the 
left-hand part of the phase plane in figure (3) and continue to collapse with different rate. 
The ones that start from above go through the right-hand part of the phase plane. In this 
case the collapse turns into a slow expansion along the expansion part of solution (ASiii) . 
In both cases the further they started up on (ASiv) the further they will go along solution 
(ASiii). In both cases the collapse or the slow expansion turn into rapid expansion as the 
trajectories are repelled from (ASiii) towards the expanding part of solution (iv). In figure 
4 we show two typical sets of evolution for a, the scale factor, and x to demonstrate each of 
these histories. One can see that the time scale on which changes take place is completely 
different but eventually both reach the point of fast expansion. In both cases the scale 
factor becomes infinite in a finite time. 

Now, we turn our attention to an open universe by setting k = —1. There is an 
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elementary exact solution with P = and x = 1 (found already by Tseytlin) for which 

a = t + to = 00 • (30) 

This solution has the same form in the Einstein frame. All the trajectories in the right- 
hand of figure (5) approach this solution marked by a circle on this figure. Tracing these 
trajectories to the past reveals that all the trajectories arc repelled from the contracting 
part of the asymptotic solution (ASiv). Here since = requires 6x-P + 6%^ + -P^ — 6 = 0, 
all the trajectories start above the exact solution (iv). Because of this, trajectories that 
occupy the left-hand of figure (5) are repelled from below by the contracting part of (ASiv). 
These trajectories are attracted by the expanding part of (ASiii), i.e these universes end 
up in rapid expansion. We see that in this case, unlike the case when k = 1, there are 
two completely different types of behavior depending on the initial conditions. Figure (6) 
shows two of the line evolutions of typical example of these two different types of trajectory. 
In the first case, the Universe ends up growing linearly while in the second case it starts 
from linear contraction and which then turns into rapid expansion. 

We turn now briefiy to the case of c 7^ 0. The analogs of Equation (17) are, for k = 0, 

P = -12x^ - 15Px - - 3c 

X = 3x'+5Px + ^'+c (31) 
(now ■ = d/dt as for the k = case). As before the phase plane properties can be explored 



from observing that 



dP 

dx 3x2 + 5Px + P^ + c ' 



The boundaries of the physical allowed regions are given by = 0, for which 

6x^ + 6Px + P^ - c = . (33) 

The solutions to Eq. (33) are therefore 

P=-3x±V3x' + c . (34) 
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It seems therefore intuitively apparent that by comparing Eq. (34) with Eq. (29) that 
c > should be qualitatively similar to /c < and the effect of c < should be qualitatively 
similar to k > 0. 

It is in principle possible to extend our method to deal with the cases k ^ 0, c ^ 
simultaneously. We do not do this here because it is too complicated to be in good taste. 

Conclusions 

Our overall conclusions are that in some stringy era, the large scale properties of the 
universe are not without some interest. Let us describe the c = 0, A; = case first. The 
contracting part of the exact solutions, the dashed line in Fig. 1, repels all the trajectories. 
So the general evolution will be away from it, into the expanding phase. As one might have 
expected by time-reversal invariance, there are solutions which are the reverse of these. 
Crudely speaking, these occupy the left-hand half of figure (1) whereas the original one 
occupy the right-hand half. These solutions start from the exact solution (iii) and are 
repelled towards the exact solution (iv) and so the universe ends up in rapid expansion. 
If one starts in the upper part of the diagram, the evolution will be towards a universe 
that reaches constant size as t — > oo, as given by the attractor shown by the dotted line, 
or analytically by equation ... If however, one starts in the lower part of the diagram, the 
expansion will always lead to hyperinflation. In our mathematical idealization, the radius 
of the universe tends toward infinity in a finite proper time. However, it could be the case 
that in a realistic theory, there is an exit regime when the effective temperature of the 
universe reduces below some critical scale, presumably the GUT scale. For k = 1, c = 0, 
(or equivalently for c < and A; = 0) then for large values of x and P, the behavior is 
qualitatively similar. However, for small x and P, there are significant differences. If one 
starts in a contracted era in the upper part of figure 3, provided one does not start too 
close to the repeller, an era of contraction will be followed by expansion, which can be very 
slow for a significant time. This is called a loitering era [7] . Eventually, the expansion will 
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accelerate, again leading to hyperinflation unless some exit from the string-dominated era 
can occur. If one starts very close to the repeller, then there is again a loitering phase which 
occurs during the contraction. Eventually however, expansion and then hyperinflation will 
set in, again subject to the caveats that we are still in the stringy era. For k = —1, c = 
(or equivalently c > and k = 0), if one starts in the upper part of the picture in a 
contracting phase, then this always turns into expansion, and at late times a{t) ~ t. This 
corresponds to the flxed point in the (x, P) plane marked by an open circle in flgure 5. 
Similarly, one can start at a point in the lower part of figure 5, and again one discovers 
that there is an attractive hyperinflating solution. 

We see therefore that during a stringy era, one can get a large variety of different 
behavior depending on the initial conditions when the stringy era sets in. Furthermore, 
depending on when the stringy era ends, one can find a number of clues to how GUT era 
physics could have started off. 

There is an addition to the cosmological interest the potential in this work to resolve 
a string theory puzzle. If string theory really does describe the physics of gravitation, 
it seems that there should be a massless scalar component to the gravitational field in 
the form of the dilaton. However, as has been emphasized by Damour and Nordvedt[8], 
the dilaton is forced by the cosomological expansion to the general relativistic limit. Our 
results are no exception to this observation as can be seen from the fact that is often 
driven to infinity, which corresponds to weak coupling of the dilaton. 
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Figure Caption 

Fig 1. A phase plane digram in the (x,-P) plane for k = and c = 0. Evolution is in the 
direction of the arrows. The exact solutions (in) and (iv) are identified by the half 
dashed half dotted lines. The dashed part is a repeller and the dotted line is an 

attr actor. 

Fig 2. Two typical solutions of Log a and x ^ function of t. 
Fig 3. The same as Fig 1. for A; = 1 and c = 0. 

Fig 4. Two typical solutions of a and x a-s function of t. The straight line corresponds to the 

points marked by the circles in Fig. 3. 
Fig 5. The same as Fig 1. for A; = —1 and c = 0. 
Fig 6. Two typical solutions of a and x cts function of t. 

a. For trajectories that start near the origin. 

b. For trajectories that start away from the origin. 
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